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The Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation is known to be 
"unstable" with respect to fluctuations in gluon virtuality, transverse momentum 
and energy requiring to go beyond the leading order BFKL. Still, these instabil- 
ities point to fruitful improvements of our deep understanding of QCD. Recent 
applications to next-leading order and to saturation problems are outlined. 



1. "Instabilities" of the BFKL Equation 
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Figure 1. The QCD evolution landscape. The BFKL evolution and its limitations, i) 
to the left: non-perturbative region; ii) to the right: DGLAP evolution; iii) to the top: 
Saturation region. 
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The Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation has al- 
ready a venerable past 1 . It appears as a key tool in many recent works 
on small-x physics (in the broad sense). It is interesting to notice that its 
limitations themselves are the seeds of interesting fields of research. Let us 
discuss limitations which can be associated with the idea of "instabilities" . 

• i) Instability towards the non-perturbative regime 

It is well known that the perturbative "gluon ladder" contribut- 
ing to the BFKL cross-section is characterized by a "cigar-shape" 
structure 2 of the transverse momenta. Hence, it is difficult to avoid 
an excursion inside the near-by non-perturbative region 

(Fig. (1), to the left). 

• ii) Instability towards the renormalization group regime 

Calculations of the next-leading BFKL kernel 3 has proven that 
the inclusion of next-leading logs gives a (too) strong correction 
to the leading log result. After a resummation motivated by the 
suppression of spurious singularities 4 ' 5 , the results show that the 
resummed NLO-BFKL kernels are very similar, e.g. "attracted" 
towards the Dokshitzer Gribov Lipatov Altarelli Parisi (DGLAP) 
evolution 6 , at least for the structure functions 

(Fig. (1), to the right). 

• iii) Instability towards the high density (saturation) regime 

The BFKL evolution implies a densification of gluons and sea 
quarks, while they keep in average the same size. It is thus natural 
to expect 7 a modification of the evolution equation by non-linear 
contributions in the gluon density. Recently, the corresponding 
theoretical framework has been settled 8 ' 9 , and is based on an ex- 
tension of the BFKL kernel acting on non-linear terms. It leads to 
a transition to the saturation regime (Fig. (1), to the top). 

The main subjects of my talk will concern contributions' 1 to point (ii), 
with a discussion of the phenomcnological relevance of (resummed) NLO- 
BFKL kernels and point (iii) , with a discussion of traveling wave solutions 
of non-linear QCD equations, as being deeply related to geometric scaling 
and the transition to saturation. 



a We shall leave thepoint (i) outside of the scope of the present conference, despite 
some recent developments related to the AdS/CFT correspondence 10 and the "BFKL 
treatment" of the 4-Supersymmetrical gauge field theory 11 . 
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2. "Instability towards DGLAP" 

The "instability" of the BFKL equation's solution w.r.t. the renormaliza- 
tion group evolution is well-known 12 . Indeed, the first correction to the 
leading- log 1/x approximation of the BFKL kernel 3 is large enough to ap- 
parently endanger the whole BFKL approach. It was soon realized that 
a large part of the problem was due to the appearance of singularities 
which contradict the renormalization group properties. Hence requiring an 
harmonization between the next leading log BFKL calculations and the 
renormalization group requirements through higher orders' resummation 
leads 4 ' 5 to a possible way out of the problem. 

Let us focus 13 on the impact of these developments on the proton struc- 
ture functions and recall the parametrization of the proton structure func- 
tions in the (LO) BFKL approximation 14 : 

f -/l(|)'"^- h "'-W'«' <» 

where Fi denotes respectively Ft, Fl, G (resp. transverse, longitudinal 
and gluon) structure functions and hi are the known perturbativc couplings 
to the photon (ha = 1 for the gluon structure function) , usually called "im- 
pact factors" 15 , xlo is the the LO-BFKL kernel, a s the (fixed) coupling 
constant and a; (7) an (unknown but factorizable 15 ) non-perturbative cou- 
pling to the proton. Mellin-transforming (1) in j — 1 = oj space, one easily 
finds 

Fi = I P~ (%Y A 7T M7) Vh) , (2) 

J 2iit \QqJ lu- t^XLoh) 

and, taking the pole contribution, one has the important relation 

w = ^-^-XLttdk^)) . (3) 

Let us try and find the equivalent relation at NLO. At (rcsummcd) 
next-to-lcading order, one can similarly write b 



b Eq.(4) is already an approximation of the (still) unknown complete (resummed) NLO- 
BFKL formula, since the photon and proton impact factors arc not yet known at NLO. 
However, one expects Eq.(4) to be a phenomcnologically valid approximation containing 
the information on the NLO kernel. 
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Figure 2. 7(0;, Q ) : T/ie "effective" F2 anomalous dimension. */(uj,Q ) has been 
evaluated from four known different parametrizations. They are all compatible in the 
range .3 < ui < 1,5 < Q 2 < 100 GeV 2 , where we restrict our analysis. 



where, by construction 

d 

— X(7, u) = Xnlo(i, w) . (5) 

The function X(7, w) appears in the solution of the Green function derived 
from the renormalization-group improved small- x Bj equation 4 , xnlo(j, 
is a resummed NLO-BFKL kernel and 

^ a s (Q 2 )= [b In (Q 2 /A^)] _1 - (6) 



c The second variable of X(-y. u>) in (5) corresponds to the choice of a reference scale 
(j, — > u) = j — 1 dictated by the treatment of the Green function fluctuations near the 
saddle-point 4 . 
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Figure 3. Test of Relation ii): LO-BFKL. 



with b = 11-2/3 N f /N c . 

At large enough Q 2 j ^qcdi one can use the saddle-point appoximation 
in 7 to evaluate (4) . Assuming that the perturbative impact factors and the 
non-perturbative function rj do not vary much d , the saddle-point condition 
reads 

Xnlo{i,u) N c a s {Q 2 ) 



V QCD 



b In (Q 2 /A2 
where 7 = j(u>, Q 2 ) is the saddle-point value. 



-XNLo(l,u) 



(7) 



d We do not take into account modifications e.g. coming from powers of 7 in the prefactors 
which may shift the saddle point 4 . We thus assume a smoothness of the structure 



function integrand around the saddle-point in agreement with the phenomenology 



13 
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Figure 4. Test of Relation ii): NLO-BFKL. 



Inserting the saddle point defined by (7) in formula (4), one obtains a 
set of conditions to be fulfilled at (resummed) NLO level as follows: 
i) The Mellin transform F 2 = Ft + Fl defines: 



d 



dln(Q 2 ) 



\nF 2 (w,Q 2 ) ~ j(oj,Q 2 ) 



ii) 7 verifies 



XNLoil) 



TT LU 



a s (Q 2 )N c ' 



(8) 



(9) 



where Xnlo is a resummed NLO-BFKL kernel. 

iii) The gluon structure function (one may also choose the obervable 
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F L ) verifies, via Mellin transform: 



In(G(w, Q 2 )) = In (P 2 {w, Q 2 )) - 7 In [h T {l) + Ml)] • (10) 



We test 13 the properties i)-iii) using NLO kernels proposed in 4 , and com- 
pared with the LO kernel condition (3). 

As an example an extraction of an "effective" F2 anomalous dimension 
i), see Fig. (2), is performed 13 using different parametrisations in a range 
of uo verifying the stability with respect to cuts on unknown (smallest) or 
irrelevant (large) XBj ■ The comparison of the property ii) to the LO BFKL 
kernel is displayed in Fig. (2) and the one with a resummed NLO-BFKL 
kernels (cf. Scheme 4 3) in Fig. (2). As is clearly seen on the figures the 
Mcllin-transform analysis disfavors the BFKL-LO kernel, while it is quali- 
tatively compatible with the resummed BFKL-NLO kernel. The remaining 
discrepancies at NLO could be attributed to finite NNLO corrections to the 
kernel or to still unknown NLO contributions to the impact factors 16 . A 
systematic study of the proposed NLO kernels is thus made possible using 
the method . 

3. "Saturation instability" 

As well-known, the BFKL evolution (even including next-leading contribu- 
tions) leads to a multiplication of partons with non- vanishing size and thus 
inevitably leads to a dense medium . This may be called the "Saturation 
Instability" of the BFKL evolution. 

The back-reaction of parton saturation on the BFKL equation has been 
originally 7 described by adding a non-linear damping term. More recently, 
the evolution equation to saturation have been theoretically derived in the 
case of scattering on a "large nucleus" , e.g. when the development of the 
parton cascade is tested by uncorrelated probes 8 ' 9 . 

In the following we will focus on the solutions of the Balitsky-Kovchegov 
(BK) equations 9 where one consider the evolution within the QCD dipolc 
Hilbert space 17 . To be specific let us consider N(Y, xqi), the dipole forward 
scattering amplitude and define 



Within suitable approximations (large N c , summation of fan diagrams, 




(11) 



"Similarly, relation iii) can be looked at using the gluon structure function parametriza- 
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Figure 5. Typical traveling wave solution. The function u(t, x) is represented for three 
different times. The wave front connecting the regions u = 1 and u = travels from the 
left to the right as t increases. That illustrates how the "strong absorption" or saturated 
phase region invades the "transparency" region. 



spatial homogeneity) and starting from the Balitsky-Kovchegov equation 9 , 
it has been shown that this quantity obeys the nonlinear evolution equation 



<9yN = ax (- 



■8l) N-aN 2 , (12) 

where a — a s N c /TT, x(l) — 2i/>(l) — V'(t) — — 7) is the characteristic 
function of the BFKL kernel 1 , L = log(fc 2 /fcQ) and ko is some fixed low 
momentum scale. It is well-known that the BFKL kernel can be expanded 
to second order around 7=5, if one sticks to the kinematical regime 8aY 3> 
L. We expect this commonly used approximation to remain valid for the 
full nonlinear equation. The latter boils down to a parabolic nonlinear 
partial derivative equation: 

Y N = a X (-d L )N-aN 2 , (13) 



with 



(14) 



The key point of our recent approach 18 is to remark that the structure of 
Eq.(13) is identical (for fixed a) to a mathematical and physical archetype of 



tions. However assumptions on the perturbative make the conclusions more qualitative 
or indicating some discrepancies to be solved at NLO. 
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non-linear evolution equation for which useful tools can be applied, namely 
the Fisher and Kolmogorov-Petrovsky-Piscounov (KPP) equation 19 for a 
function u: 

d t u(t,x) = dlu(t,x)+u(t,x)(l-u(t,x)) , (15) 

which is directly related 18 to N. The equation can be generalized to many 
physical situations, including running a. 



Fixed coupling Running coupling 




log(k 2 /Q 2 ,) log(k 2 /Q 2 ,) 



Figure 6. Evolution of the reduced front profile. Fixed coupling: left; Running coupling: 
right. The reduced front profile (k 2 /Q 2 )^ c N(k/Q S (Y),Y) is plotted against log(k 2 /Q 2 s ) 
for different rapidities. The various lines correspond to rapidities from 2 (lower curves, 
full line) up to 10 (upper curves). Note the similarity of the wave fronts, but the quicker 
time evolution (in \/t) for fixed coupling, by contrast with the slow time evolution (in 
tV3) for the running coupling case. 

Our main results are the following. The well-known geometric scaling 
property 20 is obtained for the solution of the non-linear equation (13) for 
the gluon amplitude at large energy. In our notation, the geometric scaling 
property reads 

^"^KoS 7 ))' (16) 

where Q 2 S {Y) is the saturation scale. We prove that geometric scaling is 
directly related to the existence of traveling wave solutions of the KPP 
equation 19 at large times. This means that there exists a function of one 
variable w such that 



u(t,x) ~ w(x — m(t)) 



(17) 
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uniformly in x. Such a solution is depicted on Fig. (5). The function m(i) 
depends on the initial condition. For the QCD case 18 , one has to consider 

m(t) =2t- f logf + 0(l) . (18) 

When transcribed in the appropriate physical variables, this mathemat- 
ical result, implies directly the known geometric scaling properties 18 . It is 
interesting to note how the traveling wave solutions provide a particularly 
striking mathematical realization of an "instability" as depicted in Fig. (5), 
when a stable fixed point (strong absorption) "invades" an unstable one 
(transparency). 

This mathematical analysis can be extended 18 to the study of the tran- 
sition towards geometrical scaling , i.e. the formation of the front wave as 
a function of time, both for fixed and running a, see Fig. (6). 

4. Conclusion 

In the present contribution, we have discussed some aspects of the "insta- 
bilities" of the BFKL equations, i.e.: 

• i) Instability towards the non-perturbative regime 

• ii) Instability towards the renormalization group regime 

• iii) Instability towards the high density (saturation) regime 

At first sight, these instabilities could have appeared as drawbacks of 
the whole approach. On the very contrary, we have seen that the extensions 
of BFKL equation raised up by the treatment of "instabilities" appear to 
be the building blocks of most interesting recent developments towards a 
better understanding of QCD dynamics. As an example, I chose to present 
some personal recent contributions to this discussion, which are far from 
giving an idea of the whole extent of the works f which attack the problem 
nowadays. 

As a brief outlook, let us mention: 

About Point (i), not discussed here, let us mention the formal but in- 
formative discussion on the N = 4 supersymmetric QCD field theory and 
the AdS/CFT correspondence 11 . 

Point (ii): It is the subject of a developing activity which will allow 
to master the rather high technicality of the BFKL-NLO calculations and 



I mentioned quite a few of them in the reference list but I want to apologize for the 
authors and studies which I may have forgotten in this necessarily shortened review. 
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thus to penetrate the subtle aspects of the compatibility between BFKL 
and DGLAP evolution equations. 

Point (iii): Saturation with both its phenomenological and theoretical 
aspects will certainly retain the attention of the Particle Physics commu- 
nity. The challenge here is the quest for a new phase of intense QCD fields 
and the undersatanding of its dynamical properties. 
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